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Abstract 

Metastability is observed when a physical system is close to a first order phase transi- 
tion. In this paper the metastable behavior of a two state reversible probabilistic cellular 
automaton with self-interaction is discussed. Depending on the self-interaction, com- 
peting metastable states arise and a behavior very similar to that of the three state 
Blume-Capel spin model is found. 
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1. Introduction. 

Metastable states are observed when a physical system is close to a 
first order phase transition. Well known examples are super-saturated va- 
por states and magnetic hystereses [1]. In the Figure 1 the isotherms of a 
ferromagnet are depicted on the left; T denotes the temperature, m the 
magnetization, i.e., the density of total magnetic moment, h the external 
magnetic field, m* > the spontaneous magnetization, and T c the Curie 
temperature. At temperature higher than T c the magnetization is zero for 
h = 0; it is said that the system is in the paramagnetic phase [2]. Below the 
critical temperature at h = the system can exhibit the not zero values 
m* and —m*; it is said that the system is in the ferromagnetic phase. For 
h = 0, when the temperature reaches the critical value T c the system un- 
dergoes a continuous (second order) phase transition; the name is justified 
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since the order parameter m varies continuously when T c is crossed. 

In the graph on the right in Figure 1 the behavior of the ferromagnet 
at T smaller than T c is illustrated. When h = the system jumps from the 
positive magnetization phase to the negative magnetization one, or vice- 
versa; the transition is called first order since the order parameter m, which 
is the first derivative of one of the thermodynamical potential, undergoes 
an abrupt variation [2]. Sometimes, provided the value h = is crossed 
sweetly in the experiment, the system persists in the same phase and the 
hysteresis in the picture is observed. It is then said that the phase with 
negative (resp. positive) magnetization is metastable for T < T c and h > 
(resp. h < 0) small. 




Fig. 1. Isotherms of a ferromagnet on the left; ferromagnetic hysteresis on the right. 
The temperature is denoted by T, m is the magnetization, h is the external magnetic 
field, m* is the spontaneous magnetization, T c is the Curie temperature. 

The rigorous mathematical description of this phenomenon is relatively 
recent. Not completely rigorous approaches based on equilibrium states 
have been developed in different fashions. The purely dynamical point of 
view revealed more powerful and leaded to a pretty elegant definition and 
characterization of the metastable states; the most important results in this 
respect have been summed up in [1]. 

In this paper we stick to the dynamical description and investigate com- 
peting metastable states. This problem shows up in connection with many 
physical processes, such as the crystallization of proteins [3] and in glasses, 
in which the presence of a huge number of minima of the energy landscape 
prevents the system from reaching the equilibrium [4]. The study of these 
systems is difficult, since the minima of the energy and the decay pathways 
between them change when the control parameters are varied. It is then of 
interest the study of models in which a complete control of the variations 
induced on the energy landscape by changes in the parameters is possible. 
In Section 2 we discuss the metastable behavior of the Blume-Capel model 
relying on results in [7]. In Section 3 the obtained result will be compared 
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with the known metastable behavior of reversible Probabilistic Cellular au- 
tomata with self-interaction. 

2. The Blume Capel model. 

The Blume-Capel model has been introduced in [5,6] in connection 
with the liquid Helium transition. In the context of metastability this model 
revealed very interesting for the three-fold nature of its ground states, see [7, 
8]. Consider the two-dimensional torus A = {0, . . . , L — l} 2 , with L even, 
endowed with the Euclidean metric; x,y £ A are nearest neighbors iff their 
mutual distance is equal to 1. Associate a variable a(x) = 0, ±1 with each 
site x € A and let Q = {—1, 0, be the configuration space. The energy 

associated to the configuration a € f2 is 

(1) H(a) = Y, W*) ~ <v)f ~ A 5>(*)) 2 " h E < x ) 

<x,y> xeA xeA 

where < x,y > denotes a generic pair of nearest neighbors sites in the 
torus A, A € M is the chemical potential, h € R is the external magnetic 
field, and \h\, |A| < 1. The function H will be also called Hamiltonian. 
The equilibrium behavior of the system is described by the Gibbs measure 
fx(a) := exp{— f3H(a)}/Z , where is the inverse of the temperature and 
the normalization constant Z is called partition function. 

It is possible to introduce the stochastic version of the model by defining 
a serial dynamics reversible w.r.t. the Hamiltonian (1). It will be a discrete 
time Glauber dynamics, that is a Markov chain with state space Q, and 
transition matrix p : $7 x Q, — > [0, 1] such that 

(2) p(a,r]) := J— e -P™ x {H( v )-H(a),o} 

2|A| 

for a, r] € f2 such a and r) are nearest neighboring configurations, i.e., a is 
equal to r\ excepted for the value of the spin associated to a single site; 
p(o,r\) := for a, i] € O such that a ^ i] and a and r\ are not nearest 
neighboring, that is to say they differ for the values of the spins associated 
to at least two sites. To ensure the correct normalization of the transition 
matrix, we also set pip, a) = 1 — Yl-q^a P( a > r ?) f° r an y ° e 

This dynamics, called Metropolis algorithm, satisfies the two following 
important properties: (i) only transitions between nearest neighboring con- 
figurations are allowed; (ii) the dynamics is reversible w.r.t. the Hamiltonian 
(1), i.e., 

(3) l^(o-)p(a, i]) = fJ,(rj)p(ri, a) 
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for any a, 77 € 0. The equation (3) is called detailed balance condition. 

The definition (2) of the dynamics implies that transitions decreasing 
the energy happen with finite probability," while transitions increasing the 
energy are performed with probability tending to zero for (3 — ► oo, that is 
when the temperature tends to zero. This means that when the tempera- 
ture is small, the system takes a time exponentially large in (3 to leave a 
local minimum of the Hamiltonian, i.e., a configuration a € f2 such that 
H(ri) > H(a) for any r\ € nearest neighbor of a. We can then expect that, 
wherever started, the systems tends to reach the ground state of the energy, 
i.e., the minimum of H, in a tunneling time depending on the initial condi- 
tion. Supposing that there exists initial data for which the tunneling time is 
exponentially large in j3, it is rather natural to define the metastable state 
as the configuration to which corresponds the maximum tunneling time. 




Fig. 1. Definition of metastable states. 

More precisely, following [1] and referring to the Figure 1 for a de- 
scription of the following definitions, given a sequence of configurations 
u) = uji, . . . ,u! n , with n > 2, we define the energy height along the path 
<jj as = max£ = i it,,] H(u)i). Given A, A' C 0, we let the communica- 
tion energy $(A, A') between A and A' be the minimal energy height <f> u 
over the set of paths u starting in A and ending in A'. For any <j 6 !1, we 
let I ff C fl be the set of configurations with energy strictly below H(a) 
and V a = <J>(<7, I a ) — H(cr) be the stability level of a, that is the energy 
barrier that, starting from a, must be overcome to reach the set of con- 
figurations with energy smaller than H(a); we set V a = oo if I a = 0. 
We denote by s the set of global minima of the energy (1), i.e., the col- 
lection of the ground states, and suppose that the communication energy 
r = maXo-gQ^Qs V a is strictly positive. Finally, we define the set of metastable 
states f2 m = {77 € : = T}. The set f2 m deserves its name, since in a 
rather general framework it is possible to prove (see, e.g., [9, Theorem 4.9]) 
the following: pick a £ J) m , consider the chain o~ n started at ctq = a, then 
the first hitting time Tf> = inf{t > : at € f2 s } to the ground states is a 
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random variable with mean exponentially large in (3, that is 



(4) 



/3->oo 




with Eo- the average on the trajectories started at a. In the considered 
regime, finite volume and temperature tending to zero, the description of 
metastability is then reduced to the computation of f2 s , T, and O m . 

After this rather general discussion on the definition of metastable states 
we get back to the study of the Blume-Capel model and note that rigorous 
results have already been found in [7] in the region h > A > 0. In this section 
we review those results on heuristic grounds and extend the discussion to 
the whole region h > and h > — A. 



First of all we describe the structure of the ground states of the Hamil- 
tonian. Denote by d, u and the configurations with all the spins in A 
equal respectively to —1, +1 and 0, and remark that E(u) = —L 2 (X + h), 
E(d) = —L 2 (X — h), and E(0) = 0. It is not difficult to prove that for 
A = h = the ground state is three times degenerate and the configura- 
tions minimizing the Hamiltonian are d, u and 0; for h > and h > —A, 
the ground state is u; for h < and h < A the ground state is d; for A < 
and A < h < — A the ground state is 0; for h = 0, A > the ground state is 
two times degenerate and the configurations minimizing the Hamiltonian 
are d and u; for h = A < the ground state is two times degenerate and the 
configurations minimizing the Hamiltonian are d and 0; for h = — A > the 
ground state is two times degenerate and the configurations minimizing the 
Hamiltonian are u and 0. These results are summarized in the graph in the 
left in Figure 2. Note, also, that E(0) > E(d) > E(u) for < h < A < 1, 
E(0) = E(d) > E(u) for < h = A < 1, and E(d) > £7(0) > E(u) for 
h > |A|, see the two graphs on the right in the Figure 2. 

The obvious candidates to be metastable states are the configurations 
d or 0; in particular the situation in the region h > A > looks really 
intriguing. In order to prove rigorously that one of them is the metastable 
state, one should compute T and prove that either Vd or Vq is equal to 




Fig. 2. Ground states of the Blume-Capel model. 
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r. This is a difficult task, indeed all the paths uj connecting d and to u 
should be taken into account and the related energy heights $ w computed. 
This problem has been solved rigorously in [7] in the region h > A > 
under the technical restriction 2(h/X) 2 + h/X — 1 < 2J/A. There it has 
been proven that the metastable state is d and that, depending on the 
ration h/X, during the tunneling from the metastable to the stable state 
the configuration is visited or not visited. 

As mentioned above we develop an heuristic argument to characterize 
the behavior of the system in the whole region h > and h > —A. To 
characterize the local minima of the Hamiltonian, it is necessary to compute 
the energy variation under the flip of a single spin. Then consider a € £1, 
x € A, a € {— 1,0, +1}, and denote by a% the configuration such that 
a x(y) = a (y) f° r all y 7^ x and c"(x) = a; note that a% = a iff a = <r(x). 
By using (1) we easily get 

(5) H{a a x )-H{a) = -2{a - a(x))S a (x) - (A- 4)(a 2 -a{xf) - h(a - a(x)) 

where S a {x) is the sum of the four spins of a associated to the nearest 
neighbors of the site x. Equation (5) can be used to compute the energy 
difference involved in all the possible spin flips; the results are summarized 
in the Table 1. Note that the three cases not listed in the table can be 
deduced by changing the sign accordingly, for instance if c(x) = — 1 and 
a = +1, we get H(a x ) — H(a) = — AS a {x) — 2h whose sign is positive for 
So-(x) < —1 and negative for S a (x) > 0. It is also worth remarking that 
the results on the sign of the energy differences listed in the third column 
of the Table 1 strongly depend on the assumption |A|, \h\ < 1. 



cr(x) 


a 


H{a%) - H(a) 


sign 


+ 1 


-1 


AS a {x) + 2h 


> if S a {x) > 
< if S a {x) < -1 


+ 1 





2S a {x) -4 + X + h 


> if S a (x) > +2 
< if S a (x) < +1 





-1 


2S a (x)+A-X + h 


> if S a (x) > -1 

> if S a (x) > -2 and h > X 

< if S a {x) < -2 and h < X 

< if S a (x) < -3 



From the results in Table 1 it follows that for h > X the local configura- 
tions in which a minus can appear in a local minimum are those such that 
the sum of the neighboring spins is smaller than or equal to —3, see the two 



6 



DOI: 10.1685/CSC09XXX 



configurations on the left in Figure 3. For h < A the local configurations 
in which a minus can appear in a local minimum are those such that the 
sum of the neighboring spins is smaller than or equal to —2, see the four 
configurations in the Figure 3. 

-00 + 

Fig. 3. Minus spins allowed in a local minimum; for h > A only the two configurations 
on the left are allowed, while for ft < A all the four depicted configurations are possible. 

From the first two lines in Table 1 it follows that the sole local con- 
figurations in which a plus spin can appear in a local minimum are those 
such that the sum of the neighboring spins is greater than or equal to 2, 
see Figure 4. 

+ 0-0 
+++ +++ +++ 0++ 

+ + + + 

Fig. 4. Plus spins allowed in a local minimum. 

We discuss in detail the case h > A; the analogous results in the region 
A > h > will be summarized in the Figure 6. From the necessary con- 
dition for a minus in a local minimum, see the two graphs on the left in 
the Figure 3, we have that for a configuration to be a local minimum it is 
necessary that the zeroes form well separated rectangles possibly winding 
around the torus. To verify that this condition is sufficient for the config- 
uration to be a local minimum we note that, in this case h > A, the local 
configurations in which a zero can appear in a local minimum are those 
such that the sum of the neighboring spins is greater than or equal to —2 
and smaller than or equal to +1. In the Figure 5 the possible local configu- 
ration for a zero with at least a neighboring plus are shown. This condition 
is surely met in a configuration in which the zeroes form separated rectan- 
gular clusters plunged in a sea of minuses with side lengths larger or equal 
to two. Moreover, see the Figure 3, in a local minimum direct interfaces 
between minuses and pluses are forbidden, then the pluses must necessarily 
be located in the bulk of the zero rectangular droplets. From the results in 
the Figure 5, see in particular the two graphs on the right, it follows that 
the pluses must a form well separated rectangular clusters, possibly wind- 
ing around the torus, inside a rectangular zero cluster. Note that the plus 
cluster can be separated by the minus component even by a single layer of 
zeroes. 
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- + + - - - - - 

+ + + + + + 

Fig. 5. Zero spins with at least a neighboring plus allowed in a local minimum for h > A. 

In order to study the nucleation of the stable state starting from the 
possibly metastable states and d the interesting local minima, in the 
case h > A, are the zero rectangular droplets in the see of minuses, the 
plus rectangular droplets in the sea of zeroes, and the frames made of a 
plus rectangular droplet plunged in the sea of minuses and separated by 
the minus component by a single layer of pluses (the frame). The local 
minima can be used to construct the optimal paths connecting d and to 
the ground state u. 

Consider, first, the paths from d to 0. Optimal paths can be reason- 
ably constructed via a sequence of zero droplets. The difference of energy 
between two zero droplets with side lengths respectively given by £, m > 2 
and £,m + l is equal to 2 — (h — X)£. It then follows that the energy of a such 
a droplet is increased by adding an £-long slice iff I < [2/(h — A)J + 1 = £% 
where [^J denotes the largest integer smaller than the real x. The length 
£ d is called the critical length. It is reasonable that the energy barrier Vq 
is given by the difference of energy between the smallest supercritical zero 
droplet, i.e., the square zero droplet with side length £ d , and the configu- 
ration 0; by using (1) we get that such a difference of energy is equal to 

r°=4/(h-\). 

A path from to u can be constructed with a sequence of plus droplets. 
By using (1) we get that the difference of energy between two plus droplets 
with side lengths respectively given by £, m > 2 and £, m + 1 is equal to 
2 — 2(h + X)£. It then follows that the energy of a plus droplet is increased 
by adding an ^-long slice iff £ < [2/(h + A) J + 1 = The length £% is the 
critical length for the plus droplets; the difference of energy between the 
smallest supercritical plus droplet and is equal to TJJ = 4/(/i + A). 

A path from d to u can be constructed via a sequence of frames. It is 
not difficult to prove that the difference of energy between two frames with 
internal (rectangle of pluses) side lengths respectively given by £, m > 2 
and £, m + 1 is equal to 4 — 2(h — A) — 2h£, so that the critical length for 
those frames is given by £ f d = [(2 -{h- \))/h\ + 1 and the difference of 
energy between the smallest supercritical frame and d is equal to = 
8 + 2{£ d ) 2 h - 4h£ f d e - 4(h - A), where e = £ f d -[(2-(h- \)/h]. 

Remarked that for h, A <C 1 one has rjj ~ 8/h, by comparing the energy 
barriers computed above, it is possible to find the communication energy T 



8 



DOI: 10.1685/CSC09XXX 



and to deduce all the results summarized in the Figure 6. 




rQ_ 4 r u — 8 

'd-fc-A ' d— h 



a 00 »□ 



r d >r d >r^ => o»={d} 

exit path = d — -0 — -u 
exit time = exp{ — rJj/7"} 




r0_ 4 r u_8 
1 d~ h-X 1 d~ h 



A-fi 

r d >r d , rs>r° n'»={0} 

exit path — - n 

exit time — exp{ — Tq/7 - } 



f-u 4 

1 — /,+ A 



1 d~ h-A 'd-/, 


|-u__4_ 






oQ 


u 2 

A-/7 


2-ft+A 


2 

A+/j 


rS>r d , r^>rs => n»={d} 






exit path = d — m 






exit time = exp{ — rjj/T} 








[ b: 2A>h>A>0 





00 ,0 



d|0 

2-jH-A 



rg>r^, rj>rg =► n»={d} 

H(0)-H(d)~|A| => exit path = d^u 
exit time — exp{ — fj/T} 



a A ■)■ -0 



Fig. 6. Summary of results for the Blume-Capel model. 



3. Probabilistic cellular automata with self interaction. 

We have seen above how in the case of a three-state model as the Blume- 
Capel model competing metastable states shows up. In some sense this 
result is natural because the single site configuration space is three-state. 
In the framework of Probabilistic Cellular Automata it has been shown, 
see [10,11,12,13], how competing metastable states arise in the context of 
a genuine two-state model. 

Consider the two-dimensional torus A = {0, . . . , L — l} 2 , with L even, 
endowed with the Euclidean metric. Associate a variable cr(x) = ±1 with 
each site x G A and let £1 = {— 1,+1} A be the configuration space. Let 
(3 > and k, h G [0, 1]. Consider the Markov chain a n , with n = 0, 1, . . . , 
on f2 with transition matrix 

(1) p(a, rj) = \\ p x>(7 (rj(x)) Ma, rj G Q 

where, for x G A and a G fi, p x ,a{') is the probability measure on {—1, +1} 
defined as p x ,<r{s) = !/[! + exp {— 2f5s(S a (x) + h)}] with s G {—1, +1} and 
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s a(x) = J2 y eA K ( x -y) a (v) where K(x-y) is if \x-y\ > 2, 1 if \x-y\ = 1, 
and k if \x — y\ = 0. The probability p x ,a(s) for the spin a(x) to be equal to 
s depends only on the values of the spins of a in the five site cross centered 
at x. The metastable behavior of model (1) has been studied in Ref. [11] 
for k = and in Ref. [10,12] for k = 1. 

The Markov chain (1) is a probabilistic cellular automata (PCA); the 
chain a n , with n = 0, 1, ... , updates all the spins simultaneously and in- 
dependently at any time. The chain is reversible with respect to the Gibbs 
measure fi(a) = ex.p{—(3H(a)}/Z with Z = X^en ex P{ — PH(rj)} and 

(2) H{a) = -h^2 <r(x) ~ \ E lo § cosh ^ + h )\ 

xeA 1 xeA 

that is detailed balance p(a, rf) e~^ H ^ = p(rj, a) e~~P H ^> holds and, hence, 
fi is stationary; 1/(3 is called the temperature and h the magnetic field. 

Although the dynamics is reversible w.r.t. the Gibbs measure associated 
to the Hamiltonian (2), the probability p(o~, rj) cannot be expressed in terms 
of H(a) — H(r)), as usually happens for Glauber dynamics. Given a, r\ G $7, 
we define the energy cost 

(3) A(^) = - lim = Wat) + h\ 

1 i£A: 

ri(x)[S<r(x)+h]<0 

Note that A(cj, 77) > and A(cr, rj) is not necessarily equal to A(rj,a); it 
can be proven, see [12, Section 2.6], that 

(4) e -/JA(<r,u)-/87(0) < p fa v ) < e -PA(a,r,)+f3j(P) 

with 7(/3) — > in the zero temperature limit (3 — > 00. Hence, A can be 
interpreted as the cost of the transition from a to r\ and plays the role that, 
in the context of Glauber dynamics, is played by the difference of energy. 

In this context the ground states are those configurations on which the 
Gibbs measure \x concentrates when f3 — > 00; hence, they can be defined as 
the minima of the energy 

(5) E{a)= lim H(a) = -h J^ a ( x ) 1-^0) + h\ 

^ °° x-gA x-gA 

For X C £1, we set E(X) = min a& x E(o). For h > the configuration 
u, with u(x) = +1 for x G A, is the unique ground state, indeed each site 
contributes to the energy with —h—(4 + K+h). For h = 0, the ground states 
are the configurations such that all the sites contribute to the sum (5) with 
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4 + k. Hence, for k £ (0, 1], the sole ground states are the configurations u 
and d, with d{x) = —1 for x £ A. For k = 0, the configurations c c , c° 6 
such that c e (x) = (— 1) X1 + X2 and c°(x) = (— i^i+^+i f or x _ (x\,X2) G 
A are ground states, as well. Notice that c e and c° are chessboard-like 
states with the pluses on the even and odd sub-lattices, respectively; we 
set c = {c e , c°}. Since the side length L of the torus A is even, then E(c°) = 
E(c°) = E(c). By studying those energies as a function of k and h, recalling 
that periodic boundary conditions are considered, we get E(u) = — L 2 (4 + 
k + 2h), E(d) = -L 2 (4 + K-2K), and E(c) = -L 2 (4 - «); hence E(c) > 
E(d) > E(u) for < h < k < 1, E{c) = E(d) > E(u) for < h = k < 1, 
and E(d) > E(c) > E(u) for < k < h < 1. 

In [13] the metastable behavior of this model has been studied with 
an heuristic argument very similar to the one developed in the Section 2 
to discuss the metastable behavior of the Blume-Capel model. For the 
details we refer the interested reader to the quoted paper, we just mention 
here, that quite surprisingly results very similar to the ones obtained in 
the framework of the Blume-Capel model are found, provided the different 
parameters are interpreted according to the correspondences in Table 1. 



Blume-Capel 


u 


d 





h 


A 


PCA 


u 


d 


c 


h/2 


k/2 



Notice that the role of the zero state of the Blume-Capel model is 
played, in the context of the PCA, by the flip-flopping chessboard-like 
configurations. As (4) shows, the discussed results are valid in the limit 
[3 — > oo. Their validity at finite temperature can be tested with Monte 
Carlo simulations, see the configurations in Figure 1 observed in a run 
of the dynamics of the PCA with the parameters specified in the caption 
and with starting configuration d. On the left it is shown that if the self- 
interaction is present the nucleation of the plus phase is achieved directly; 
the plot on the right shows that, if the self-interaction is zero, than the 
chessboard-like phase is visited before the plus phase is nucleated. 

REFERENCES 

1. E. Olivieri and M. E. Vares, Large deviations and metastability. Cam- 
bridge University Press, UK, 2004. 

2. K. Huang, Statistical Mechanics. Wiley, 1987. 

3. P.R. ten Wolde and D. Frenkel, Enhancement of Protein Crystal Nucle- 
ation by Critical Density Fluctuations, Science, 277 (1997) pp. 1975- 



11 



E.N.M. Cirillo, F.R. Nardi, C. Spitoni 




Fig. 1. On the left, typical configuration of the PC A with k = 1; simulation performed 
on a 380 x 230 torus with j3 = 0.7. White and black points represent respectively minus 
and plus spins. On the right, the same with k = 0; the chessboard region looks grey. 

1978. 

4. G. Biroli and J. Kurchan, Metastable states in glassy systems, Phys. 
Rev. E, 64 (2001) pp. 016101(1-14). 

5. M. Blume, Theory of the First-Order Magnetic Phase Change in UO2, 
Phys. Rev., 141 (1966) pp. 517-524. 

6. H.W. Capel, On the possibility of first-order phase transitions in Ising 
systems of triplet ions with zero-field splitting Physica, 32, (1966) 
pp. 966-988. 

7. E.N.M. Cirillo and E. Olivieri, Metastability and nucleation for the 
Blume-Capel model. Different mechanisms of transition, J. Stat. Phys., 
83 (1996) pp. 473-554. 

8. T. Fiig, B.M. Gorman, P.A. Rikvold, and M.A. Novotny, Numerical 
transfer-matrix study of a model with competing metastable states, 
Phys. Rev. E, 50 (1994) pp. 1930-1947. 

9. F. Manzo, F.R. Nardi, E. Olivieri, and E. Scoppola, On the essential 
features of metastability: tunnelling time and critical configurations, J. 
Stat. Phys., 115 (2004) pp. 591-642. 

10. S. Bigelis, E.N.M. Cirillo, J.L. Lebowitz, and E.R. Speer, Critical 
droplets in metastable probabilistic cellular automata, Phys. Rev. E, 
59 (1999) pp. 3935-3941. 

11. E. N. M. Cirillo and F. R. Nardi, Metastability for the Ising model with 
a parallel dynamics, J. Stat. Phys., 110 (2003) pp. 183-217. 

12. E.N.M. Cirillo, F.R. Nardi and C. Spitoni, Metastability for a reversible 
probabilistic cellular automata with self-interaction. J. Stat. Phys., 132 
(2008) pp. 431-471. 

13. E.N.M. Cirillo, F.R. Nardi and C. Spitoni, Competitive nucleation in 
reversible Probabilistic Cellular Automata with self-interaction. Phys. 
Rev. E, 78 (2008) pp. 040601(1-4). 



12 



